Ensemble based quantum metrology 
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The field of quantum metrology promises measurement devices that are fundamentally superior to conven- 
tional technologies. Specifically, when quantum entanglement is harnessed the precision achieved is supposed 
to scale more favourably with the resources employed, such as system size and the time required. Here we 
consider measurement of magnetic field strength using an ensemble of spins, and we identify a third essential 
resource: the initial system polarisation, i.e. the low entropy of the original state. We find that performance 
depends crucially on the form of decoherence present; for a plausible dephasing model, we describe a quantum 
strategy which can indeed beat the standard quantum limit. 



Quantum metrology deals with the physical limits to mea- 
surement. Typically one prepares a probe system of size K in 
a suitable initial state, this system acquires information about 
the quantity of interest, and then the probe system is mea- 
sured. The process may be repeated either with a series of 
probes over time or, as in the present analysis, with multiple 
probe systems simultaneously. The way in which the probe 
is prepared is closely related to the uncertainty of the param- 
eter estimation: if prepared in a non-entangled state then the 
minimal uncertainty achievable scales with \j\[K [ 1 ], the so- 
called Standard Quantum Limit (SQL). However, this limit 
can be beaten if we allow arbitrary states for the prepara- 
tion of the probe system, i.e. if we include entangled states, 
as demonstrated in recent experiments J2J-2]. In idealised 
cases, the minimal uncertainty achievable scales with 1 jK - 
the Heisenberg limit [6| - which can be achieved by making 
use of GHZ-states, also called NOON-states Q. 

In optical quantum metrology the probe system is a partic- 
ular state of A' photons, for example a NOON-state, which is a 
superposition of 'all K photons in channel A' with 'all photons 
in channel B' . When an optical element inducing an unknown 
phase shift is placed in channel A, then the probe acquires an 
internal phase K times as great as that which would be ac- 
quired by a single photon, and information about this phase is 
measured through an interference effect. In the present analy- 
sis we consider an analogous experiment involving K atomic 
spins in a large molecule, which probes the strength of an ex- 
ternal magnetic field. An essential difference is that we con- 
sider a large ensemble of probe molecules which are neces- 
sarily prepared, exposed to the field, and ultimately measured 
collectively (8| [9| - that is, addressing of individual probe 
molecules is impossible. 

The dynamics of an ensemble is typically observed by mea- 
suring the free induction decay (FID) spin signal. Monitoring 
the FID can be seen as a continuous and simultaneous mea- 
surement of two non-commuting observables iflOl . Here the 
observed system is barely altered by the measurement, as the 
number of spins in a typical sample is usually so large. Given 
this type of measurement, no scaling laws for the precision 
have been reported yet. Recent studies have rather looked at 



the effect of temperature on the Fisher information of three 
types of states ifTTIl and waveform estimation and its implica- 
tions for quantum sensing [12]. 

In this Letter we compare two strategies for measuring a 
small shift 8B of a probe magnetic field from a reference field 
in a spin ensemble setting [8 9|. This problem is equivalent 



to measuring the Larmor frequency 8 



SB 



of a precessing 



spin in the probe field, where y is the gyromagnetic ratio. In 
the quantum strategy, we consider a macroscopic ensemble 
of N sensor molecules each consisting of K spins where each 
molecule is prepared in a very sensitive entangled state (see 
Fig. fl). This state senses the field for a wait time T w by ac- 
quiring a phase K8T W . We then map the phase onto one spin 
of the molecule from which it can be read-out by observing 
the free induction decay (FID). As a performance benchmark, 
we compare this to the classical or standard strategy, where 
we determine the Larmor frequency by observing the FID of 
the same number of uncoupled spins. We will focus on the 
resources consumed by the two strategies and on the impact 
of decoherence. 
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FIG. 1. (a) Schematic of a sensor molecule with five satellites, (b) 
Quantum circuit employed in the quantum strategy. 

The role of decoherence is very non-trivial. In the absence 
of decoherence and for a perfect projective measurement, the 
lowest achievable uncertainty in estimating 8 is 1 /K>/N for 
the quantum strategy and \/\/KN for the classical strategy 
1 6 1. Hence, the discrepancy in the precision of the two strate- 
gies increases with the size of the entangled state K. How- 
ever, the spins will be subject to decoherence in any real- 
world experiment, and the sensitive entangled states will de- 
cohere faster than separable states when K becomes larger 
JT3J. This increased effective decoherence rate competes with 



the enhanced precision, making it difficult to predict the per- 
formance of the quantum strategy. 

In the following, we shall determine lower bounds for the 
uncertainty of the parameter estimation from a measured FID 
for both strategies. First, we must generate a suitable en- 
tangled (GHZ) state to obtain enhanced sensitivity with the 
quantum strategy, meaning we require some amount of quan- 
tum control over the molecules. The details of how this is 
accomplished are unimportant and we shall in following con- 
sider the example shown in Fig. [Tl molecules with this star 
topology have been employed in recent experiments J8] |9) . 
Each molecule consists of one central spin A and K — 1 non- 
interacting satellites of type B. The satellites interact with the 
central spin through an Ising type interaction, leading to the 
following Hamiltonian in an external magnetic field (in natu- 
ral units, i.e. h = 1): 

k-\ K-\ 

H K = y8at + 8 £ a J " +J £ ff- A <8> a,' , (1) 



where a, denotes the usual Pauli matrix with eigenvalues ±A. 
8 denotes the Zeeman splitting of the B spins and y — 2& the 
ratio of the gyromagnetic numbers of A and B, and we assume 
7 « 1 . J describes the Ising interaction strength. 8 is the 
parameter to be estimated. 

We focus on the fundamental comparison between classical 
and quantum strategy given a fully polarized initial state of the 
sample: 

Po = |0)a(0|a<8>|0...0) b (0...0| b = |00)(00| , (2) 

where the underscore denotes the state of the K — 1 satellite 
spin register. For the example molecule shown, the GHZ- 
state can be constructed from this initial state by applying the 
pulse sequence shown in Fig. [Tl i.e. first a Hadamard gate on 
the central spin and then a CNOT gate (control qubit = cen- 
tral spin, target qubits = satellites). The resulting GHZ-state 
4= (| 00) + 111}) freely evolves for a time T w , acquiring phase 
K times faster than a single spin. However, at the same time 
it is also more vulnerable to decoherence. In practice, the de- 
phasing rate of an individual spin, a — — J» < 0, is limited 

'2 

by inhomogeneous broadening [ 14 1, thus we can neglect spin 

flip processes. Of course the dephasing rate of the GHZ-state 
j3 < is related to a, and we shall discuss this dependence 
later. The state of the system after the wait time T w is then: 



P3 = ^(|00)(00| + |ll)(ll|+ e 



* 5r "<+^|00)<ll 
K5T » i+ P T »\U)(00\ 



(3) 



To measure the acquired phase K8T W we map the GHZ-state 
onto the central spins by applying a CNOT gate. These spins 
are then measured by observing the decay of the transverse 
magnetization at M £ N discrete points in time, separated by 
the sampling time t s . The observed signal x m at time mt s (m = 



0, . . . ,M— 1) can be modelled as a sum of Gaussian distributed 
noise b m and the ideal signal x m — (X + iY) (mt s ) [ 15 ] 

X m = X m + b m = ce K8T w i+pT We S,nt s+ a,« ts + ^ ^ (4) 

where c is a proportionality factor that depends on the number 
of molecules in the sample. A simulation of the FID is shown 
in Fig. [2] 
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FIG. 2. (color online) Simulation of the free induction decay (FID) 
of the classical (upper) and quantum strategy (lower); see Eqs. ({9}) 
and (l|4j). The simulated measurement points (crosses) were cho- 
sen randomly from a Gaussian distribution with a standard deviation 
of 0.05. In the classical strategy we observe uncoupled precessing 
spins and the strength of the probe field is given by the oscillation 
frequency. In the quantum strategy, a GHZ-state senses the field for 
a time T w without producing a signal, but acquiring phase K times 
faster than a single spin and dephasing at a rate j8. The acquired 
phase is mapped onto the central spin and the FID is measured. The 
strength of the magnetic field can be now estimated from the phase 
and the oscillation frequency. 

A suitable metric for the precision of the measurement is 
the Cramer-Rao bound (CRB) [15|, which essentially offers a 
lower bound on the uncertainty (standard deviation) a„, of an 
estimated parameter pc (£= 1,...,P): 



Gpf ^ CRB P( 



:f- 1 ) 



(5) 



Here, F denotes the Fisher information 1 16 1 given by the real 
part of a complex-valued matrix product 



F = \3i(D f D) 



(6) 



where D, 



dp. 



and the partial derivatives are evaluated at 
the parameters that we are going to estimate. In our case the 
parameters are c, a, and 8. a = o r = cr, denotes the standard 
deviation of the real/imaginary part of the noise. Inverting the 
3x3 Fisher information matrix we obtain: 
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£ (KT W + mt s ) exp(2amr,) 
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(7) 



Next we determine the CRB for the classical strategy. Here 
we are given N K identical and uncoupled spins. We obtain 
the relevant signal by rotating the initial state po into the xy- 
plane with a global Hadamard gate followed by measuring 
the transverse magnetization. The resulting density matrix 
evolves as 



P4(0 = 



1 



1 

,-i8t+at 



,iSt+at 

1 



giving rise to a measured signal of the following form 



xL=c'e Smt ° i+amt °+bL 



(8) 



(9) 



Analogously to above we obtain the CRB for this signal: 

1 a' 



CRB 



5,STD — 



/M'-l 
I (mt s ) 2 exp(2amt s ) 

m=0 



(10) 



This expression is a lower bound on the uncertainty for in- 
dependent spins against which we shall benchmark the quan- 
tum strategy. Hence it may be considered as the SQL in this 
setting of spin ensemble measurements. For a high signal-to- 
noise ratio (SNR) achieved by a sufficiently large number of 
molecules N, there exists an efficient estimator which matches 
the accuracy predicted by the CRB ifTTl . Hence, we can di- 
rectly use Equations (|7]i and ( fTO} to compare the two strate- 
gies. 

First, we need to specify a fair comparison with the same 
resource allocation for both strategies. For conventional quan- 
tum metrology with projective measurements the challenging 
question of a fair resource comparison has recently been ad- 
dressed in Ref. [ 18 1. In the present case of ensemble quantum 
metrology, we shall at first allow both strategies to observe the 
full FID while consuming the von Neumann entropy of 1 • N 
spins. As we have defined it, the quantum strategy consumes 
one spin per molecule (the central spin A) in the measurement 
and all other spins remain pure. This implies that we must also 
only measure 1 -N spins instead of K- N spins for the classi- 
cal strategy, meaning that the SNR of both measurements are 
equal £ = %. 

At first sight this way of counting resources may look bi- 
ased towards the quantum case, as we do not seem to take the 
satellite spins into account. Nonetheless, our comparison is 
fair: the K— 1 satellite spins act as an antenna to pick up phase 
more rapidly, yet they are not 'consumed' (this is a direct con- 
sequence of the dephasing model of decoherence). After the 
quantum strategy is complete, the central spin is measured (its 
polarisation is lost) but all satellite spins are back in the pure 
state |0_) and could be recycled to obtain a further parameter 
estimate. The accuracy of such an estimate, made using any 
sensible protocol on these remaining (K — l)N spins, can be 
no worse than that obtained using the standard strategy. This 
observation validates the classical resource count stated in the 
previous paragraph. 



The ratio of Eq. ( 10 1 and Eq. (jTli can be approximated by 
an integral expression (assuming the sampling rate resolves 
the decay, i.e. t s < 7|): 



Roo.= 



CRB 5 , STD J"(KT w +t) 2 exp(2at)dt 



(11) 



CRB 5 , GHZ V e- 2 P T « J^t 2 exp(2at)dt 
= e^'V 1 ~ 2KaT w {\ - KaT w ) . (12) 

Whenever R^ > 1 the quantum strategy outperforms the clas- 
sical strategy with respect to the precision of the parameter 
estimation. 

So far we have not yet specified how the dephasing rate j8 of 
the GHZ-state relates to the dephasing rate a of a single spin, 
and we shall now discuss a number of different decoherence 
models for the GHZ-state. 

First, we consider j5 = a, implying that the GHZ-state 
does not decohere faster than a single spin; this is expected 
for a given macroscopic magnetic field inhomogeneity. In 
this case, there is for any K > 1 a wait time T w for which 
the quantum strategy surpasses the classical one. Conversely, 
completely correlated or collective noise over each molecule 
has the most aggressive effect on the GHZ-state |19|. Here 
the noise can be described with a single Lindblad operator 
Y,j=\ C z and j3 = K 2 a, while for uncorrected noise j3 = Ka 
ll20l . In general we consider a power law dependence of j3 on 
a, i.e. j3 = K p a, where < p < 2. In recent experiments the 
decoherence rates for highly correlated solid-state spin states 
were obtained experimentally [21J, revealing p sa 1/2. The 
authors attributed this to non-Markovian correlated noise. A 
significantly smaller value for p was found in [9|, where the 
T 2 * time of a single spin was determined to be 0.37 s and that 
of a 13 particle GHZ-state to be 0.28 s, which can be inter- 
preted as a factor of p ss 0. 1 1 . 




FIG. 3. (color online) The maximal £%oo for a given A" as a function 
of K and p. £%„ > 1 implies that the quantum strategy outperforms 
the standard strategy. 

The performance of the quantum sensor depends critically 
on the value of p. One can easily check that for 1 < p < 
2, R^ < 1 for any T w > and Roo = 1 only for T w = 0. We 



therefore conclude that the precision of the quantum method 
never outperforms the precision of the standard method, if 1 < 
p < 2. In contrast, f or p < 1 there is an optimal wait time 
for T w the GHZ state which is non-zero; this varies with K 
but always gives R oa > I. Specifically, calculating the optimal 
waiting time T w that gives rise to a maximum of ^ m (for a 
given K and p < 1) yields 



maxRo, 



exp (\ (\ -KP- l +Vl -K 2 P~ : 



(13) 

which scales like ^K p ~ l to leading order (see Fig. 31. There- 
fore the ensemble SQL can indeed be beaten with a quantum 
strategy if the decoherence of the GHZ-state is not too aggres- 
sive, i.e. for p < 1 . Moreover, we see that under this condition, 
the precision of the estimation improves monotonically as K 
increases. 

So far we have neglected time as a resource, focusing in- 
stead on system size and the consumption of initial polarisa- 
tion (maximising the von Neumann entropy of one spin per 
molecule). It is interesting to extend our analysis to a re- 
stricted process time, meaning only a part of the FID can be 
observed. Since the first part of the FID contains most infor- 
mation, this would enable a 'better' sensor by a series of re- 
peated runs in a given time window r tot if one had the ability 
to reset the spins to their initial state after time T, for exam- 
ple with an optical switch. We assume that this can be done 
instantly and then the uncertainty of the parameter estimation 
of this series is given by 



\/P 



CRB 



STD/GHZ 



= SsTD/GHZ 



(14) 



Obviously one would choose the length of a time slice T op- 
timally, i.e. in such a way that the uncertainty per a/Hz, i.e. 
Sstd/ghz IS minimal. For the classical strategy we find by us- 
ing the integral approximation for the CRB from above and 
numerical optimization that 



S* STD := minSsTD ; 



3.21y/t s \a\ 



3/2 ' 



(15) 



which is attained for the optimal time T* w 1.69r 2 *. In the 
quantum strategy 5ghz a ls° depends on T w , the amount of 
time for which the spins are in the GHZ-state. In contrast to 
the classical strategy, the minimum here depends on the sys- 
tem parameters K and p. We have not been able to find an an- 
alytic expression for Sq HZ :— min^ j w Sghz an( l therefore we 
performed a numerical optimization with the results displayed 
in Fig. |4] As in our first comparison we assume that the SNR 
for both strategies are equal. If aggressive noise is affecting 
the GHZ-state the optimal quantum strategy is basically the 
optimal standard strategy, as T* —¥ 1.69r 2 * and T* — > 0, when 
p — > 2. For small p and large K however the quantum strat- 
egy significantly outperforms the standard strategy. Interest- 
ingly the quantum strategy can now beat the optimized stan- 
dard strategy for values p that are slightly larger than one. 




FIG. 4. (color online) (a) Minimal (optimized over T and T w ) uncer- 
tainty of the parameter estimation for the quantum strategy per \/Hz, 
i.e. Sqhz m un i IS °f V^s\ a \ c m dependence of K and p. (b) 
The corresponding optimal times T* (upper surface) and T* (lower 
surface) in units of T 2 * for which the minimum in (a) is attained in 
dependence of K and p. 



In conclusion we have presented a framework for analysing 
the performance of quantum metrology using spin ensembles. 
This framework incorporates the special nature of the non- 
projective measurement process, and leads one to consider the 
polarisation (i.e. the low entropy) of the initial state as a re- 
source. We find that the decoherence model plays a defining 
role in this framework, and we have identified the parame- 
ter regime where a certain quantum strategy can beat the best 
classical strategy. 
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